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Abstract 

, For a noninteracting Bose gas with a fixed one-body Hamiltonian independent of tfie number 

^>f-^ ' of particles we derive the inequalities {Ni)N < {Ni)N+i, {NiNj)N < {Ni)N{Nj)N for i 7^ j, 

I d{N(}) N /d(3 > and (Nt)^ < {Ni)N- Here Ni is the occupation number of the ith eigenstate of 

. H^, (3 is the inverse temperature and the superscript + refers to adding an extra level to those of 
0^ ' ■ The results follow from the convexity of the A'^-particle free energy as a function of A'^. 
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PACS numbers: 05.30.Jp, 02.50.Cw 



^ |i Correlation inequalities have been playing a particular role in statistical physics. They have 

always been limited to special models, most of them being valid for classical spins on lattices 
12] .Nevertheless, when they could be applied, they proved to be very useful in deriving rigorous results, 
. ^ for example, on phase transitions or on the existence of limits of correlation functions, 

y . In this paper I present correlation inequalities for a family of continuous or lattice quantum systems: 

^ I noninteracting Bose gases in boxes, on tori or in some sufficiently fast increasing confining external 

potential. Such a system is characterized by a one-body Hamiltonian having a discrete spectrum 
£0 < £1 <i £2 < • • •■ In the continuum the number of eigenstates is infinite and we have to impose the 
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X .. 

^ . condition that the sum of the Boltzmann factors Xi — exp(— /Je^) is finite. 



tre-''^°=^e-''^'^^a;, <oo. (1) 

i=0 4=0 

The relevant random variables are the occupation numbers Ni of the eigenstates of . In the 
canonical ensemble they are weakly dependent, their only inter-dependence coming from the constraint 
Ni = N . Hence, their joint probability distribution is 

3>0 {nj}:^n,=Nj>0 

Note that Zq ^ 1. Let be a one-body Hamiltonian whose spectrum specTJ^j!. = specH^ U {e}. 
Below {■)n and (•)^ denote mean values in the A^-particle canonical ensembles generated by and 
respectively. 

Theorem, (i) {N,)n < {N,)n+i 

(ii) d{N,)Nld\nx, = {N,N,)n - {N,)n{N,)n < for i ^ j 

(iii) d{No)N/dp > 

(iv) {N,)+ < {N,)n. 
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The content of these inequahtics is intuitively obvious but not more obvious than that of the 
inequahties of Griffiths, Kelly and Sherman (GKS) which say that in ferromagnetic Ising models the 
spins are positively correlated [310]. Having in mind the subtleties of the proof of the GKS inequalities, 
one cannot expect a genuinly simple proof for this theorem either. 

Inequality (iv) can be put in a more general form. Let be a one-particle Hamiltonian, suppose 
that exp{—f3H^) is trace class and speci?° C specH^ where repeated eigenvalues are considered sep- 
arately. If Si is a common eigenvalue then the occupation number Ni of the corresponding eigenstate 
(which may be different for i?" and H^) satisfies 



where denotes the canonical expectation value with respect to H^. One obtains (iv') by a repeated 
application of (iv) or directly by a simple modification of the proof of (iv). One may wonder about the 
relevance of the fourth inequality. One-particle Hamiltonians with a modified spectrum can serve as 
auxiliary tools, as for instance in OE). Actually, the motivation of the present paper was the need of 
this inequality in proving the occurrence of a generalized Bose-Einstein condensation in an interacting 
trapped Bose gas, when both the interaction and the one-dimensional harmonic trap potential are 
scaled as N tends to infinity 

The basic ingredient of the proof of the theorem is the convexity of the iV-particle free energy. 

Lemma. The N -particle free energy is a convex function of N , namely 



Introducing Fj^ = — InZyv, © is equivalent to F^+i ~ F„i < Fn+i — Fn which is just convexity. 
The usual statement about the convexity of the free energy in homogenous systems is that in the 
thermodynamic limit the free energy density is a convex function of the particle number density. 
This holds quite generally true, whenever there is asymptotic equivalence between canonical and 
grand-canonical ensembles. Indeed, in this case the free energy density is the Legendre transform 
of the pressure which is trivially convex as a function of the chemical potential. Inequality ((SJ is 
more model dependent. It perhaps remains true in the so-called diagonal model of a Bose gas [J] if 
the pair interaction has a nonnegative Fourier transform, because then there is a repulsion between 
the occupation numbers of different plane-wave states. Also, in classical lattice models with repulsive 
interactions JSJ may hold true because dividing a finite piece of a lattice into two parts, the more equal 
the sizes of the parts the larger the number of repulsive links to be cut (and the product of the partition 
functions of the two parts with it). Note that the subadditive property InZm+n < ln.^n + ^^Zm is a 
special case of © {Zq = 1). 

Proof of the theorem. 

We prove the first inequality by using the lemma. The scheme of the proof of the rest will then be 
(i) =J> (ii) =^ {Hi) and (i) => (iv). 

(i) Let us introduce the notation Z^^i for the iV-particle partition function in a system where the 
level i is missing. 



< Zm+iZn for any m < n. 



(3) 




(4) 



The lemma applies to these partition functions as well. For any i 



N 




(5) 



Using lO and the corresponding identity for Zn+i, 

, N+l 



1 



N 





l — m.i 



m—1 



'N — 171,1 
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AT AT 



fc=0 m=l 

Af JV 



>fc=0 m=l 

AT AT 



+ 15^ — 7^ { / , 'm-Xj^[ZN,iZN+l-m,i — Zn -m.,iZN+l,i] 

+{N + l)xf+^ZN,i +Y,{N+l- m)x'^+''+^ZN-m,i \ . 



The expression in the last Hne is strictly positive and. due to the lemma, in the penultimate line each 
term of the single sum is nonnegative. Let us rewrite the double sum. Since the diagonal terms vanish 
and the difference in the square bracket changes sign if k and m are interchanged, we obtain 

X^~^"^{m — k)[ZM-k,iZN+l-m,i — ZM-m,iZN+l-k,i\- 

l<k<m<N 

Because the smallest among the four subscripts is N — m, we find again with the lemma that each 
term is nonnegative. This proves inequality (i). 
(ii) For i ^ j let us introduce 

1 " 

{Ni)n,3 = 51 kxlZn_k,i,j , (6) 

fc=l 

the mean value of Ni in the n-particle canonical ensemble of a system in which the level j is missing. 
Here Zn-k,i,j is the (n — fc)-particle partition function in a system with missing levels For any 
N > n and m = N — n, {Ni)n.j agrees with the conditional expectation value of Ni in the A^-particle 
system provided that Nj = m. The inequality (i) is valid also in this case. One can write 

N 

{Ni)N = ^Pm{Xj){Ni)N-m,j (7) 
m=0 

where 

Pm{Xj) = (8) 

are probabilities, J2m=oP"^i^j) = ^ (^"^^ Y.Z=i'^P"^i^j) = (^i)Af)- Now 



dpmjXj) _ xfZN-m,j 



9 In a; 



■j Z.JV 



(m - {Nj)N) = Pm{xj){m - {Nj)N) (9) 



changes sign when m passes through {Nj)^- Therefore, by applying the inequality (i) to {Ni)N-rn,j 
we find 



m=0 



3 



0<m<(A'j>jv ^ {Nj)N<m<N ■' 



0<m<{Nj)N ■' (Nj)N<m<N ■> 

d ^ 

< {N.)N^UN,).Ug^YP^r.ix,)=0 (10) 

^ m— 

which proves the second inequaUty. The identity 

N 

{NiNj)N = Y n^PrniXj){Ni)N-m,j (H) 
m— 1 

can be read off from the first line of H1(J|1 . 

(iii) 



oo 

= - E(^J- - eo)[(A^oiVj)Ar - {No)n{Nj)n] > (12) 

where we used 

oo 

Y[{NoNj)n - {No)n{Nj)n] = N[{No)n - {No)n] = (13) 

in the last equahty and (ii) in the inequaHty. 

(iv) Let X = exp(— /3e) and let denote the iV-particle canonical partition function of the system 
with the additional level. In analogy with Q and by applying (i) 

= E Jl'^^ mN-n. ^ E Plnix){N.}N-n. < (N,) ^ (14) 

m=0 N m=0 

because Pq{x) = Z^/Z^ < 1. 

This finishes the proof of the theorem. 

Proof of the lemma. 

Let pi < p2 < ■ ■ ■ < pi he a, partition of N, that is, pi > I and Yl\ pi — N. We denote this by 
p — {pi,p2, ■ ■ ■ ,pi) H N. For each m G [0, A^], ZmZN-m is a sum over the same set of monomials 
of xi, X2, ■ ■ ., namely all those Y[j Xj^ with rij = N. For a given m the coefficient of a monomial 
depends only on the set of nonzero nj and not on their order or on the set of the corresponding 
subscripts j. Therefore 

E«(piH E n^? (15) 

p\-N {"j}~P 

where {rij} ~ p means that, after rearranging in increasing order, the vector of nonzero Uj agrees 
with p. We shall prove that for any 1 < Pi < ■ ■ ■ < pi 

a{pi,p2, ■ ■ ■ ,pi\m) < a{pi,p2, ■ ■ ■ ,Pi\m+ 1) if m < ^{pi -\ hp;) (16) 
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which impUes ||3J) through a term- by-term inequality. The proof is obtained by induction over I. 

a(pi, . . . ,pz|m) is the number of different mth degree monomials within x^^ ' ' '^T more 
popular terms, the number of possibilities to distribute m Euros among I people in such a way that 
they can obtain at most pi, . . . ,pi, respectively. Thus, a{p\m) = 1 for any p > 1 and any < m < p, 
so that for I = 1 Hl()|) holds with equality. The induction hypothesis is that for an / > 1 

a(pi, . . . ,p/_i|to) < a(pi, . . . ,p;_i|m + 1) if m < i(pi H \-pi-i)- (17) 

Because of the identity 

a{p\m) ^ a{p\^Pi-m) (18) 

the hypothesis (|17|l is equivalent to 

a(pi, . . . ,p/_i|to) < a(pi, . . . ,p;_i|n) if m < min{n, pi H hp;-i-n}. (19) 

Let us introduce the notations 

Qipi, ...,pi) = {qeZ' ■.0<q,<p,, l<i<l} , (20) 

C{pi, . . . ,pi\m) = {q e Q{pi, . . . ,pi) -.^qi^m} (21) 

1 

i-i 

PC{pi,. . . ,pi\m) = {qe Q{pi,...,pi-i) ■.m-pi<^qi<m}. (22) 



and 



We have 



a(pi, ...,pi\m) = |C(pi, . . . ,_pz|to)| = |PC(pi, . . . ,p;|m)| (23) 



1^1 meaning the number of elements of the set A. The first equality is obvious. The second comes from 
the fact that PC{p\m) is the orthogonal projection of C(p\m) onto the I — 1 dimensional subspace 
perpendicular to the Ith unit vector and, because of the constraint qi = m, the projection is a 
bijection. Comparing (|22|) with 



i-i 



PC{pi,...,pi\m + l) = {qe g(pi,...,p,_i) : TO + 1 - p; < < m-f 1} (24) 

1 

we find that 



i-i 



\PC{p\m + l)\PCip\m)\ = |{qeQ(pi,...,pz_i) :^g,;=m+l}| 



1 



a{pi, . . . ,pi^i\m+ 1) (25) 



while 



\PCip\m)\PC{p\m+l)\ = \{qeQ{pi,...,pi^i):Y,<l^=m-Pi}\ 

1 

= a(pi, . . . ,pi_i|m -p;). (26) 

Note that either of H25|l and H26() can vanish if to -I- 1 > pi + ■ ■ ■ + pi-i or m ~ pi < 0, respectively. 
The inequality (|16|l is satisfied if 

a{pi, . . . ,pi-i\m-pi) < a{pi, . . . ,p;_i|to + 1) , (27) 
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and this follows from the induction hypothesis if 

m ~ pi < min{TO +1, pi + ■ ■ ■ + pi-i — (m + 1)}. (28) 

If 771+ 1 < (pi + - • •+p/_i)/2 then this is the minimum on the right-hand side of 128(1 which is therefore 
fulfilled. If 777 + 1 > {pi + • • ■ + p;_i)/2, the minimum is pi + • • ■ + pi-i — {m + 1). But then 

m- pi < pi H \-pi-i - (777 + 1) 

is equivalent to 

11/ 

which holds true because 77i < {pi + ■ ■ ■ +pi)/2. By this we have finished the proof. 

Concerning the proof of the lemma, two remarks are in order. First, it is tempting to present 
a more elegant argument which refers to the convexity of the hyper-rectangle Q{p)- If we had 
instead of in the definition of Q(p), C(p|777) would be the intersection of the rectangle with the 
hyperplane ^ qi = m. Because of the convexity of Q{p) the I — 1 dimensional Lebesgue measure of 
C(p|777) cannot have a minimum in the interval < 777 < ^p^. By symmetry, it has a not necessarily 
strict maximum at m = '^pi/2, implying that it is nondecreasing for m < ^pi/2. However, we 
do not need this result for the Lebesgue measure of the cut but for the number of points of integer 
coordinates on it, which makes the convexity argument somewhat shaky. 

Second, the number in question can be given by a formula. It reads 

Pi min{P3,m-Ei=l 9i} minfPi- 1 9i } 

a(p|777)=£... E 1- (29) 

91=0 ijj=inax{0, JZj <ji-5^'^^._l_^ Pi} qi_i=ma.^{0,m-Y,\z\qi-pi} 

If 777 < Pi, in all the summations the lower bounds equal zero and thus for m < pi, a(j)\m) < a(j3\m+l) 
can explicitly be seen. For m < pi by simple combinatorial considerations we find 



a(pH=E-- E — ^.-^^-'-'). (30) 

91=0 9i-i=0 



If 777 > pi the general formula is less transparent. This is why we have opted for the inductive proof. 

Acknowledgment This work has benefited from the support of the Hungarian Scientific Research 
Fund through Grant T 42914. 

Note added. After having submitted the paper, I learned from Valentin Zagrebnov that the convexity 
of the free energy was earlier shown by Lewis, Zagrebnov and Pule 8^- Their proof is somewhat 
shorter, the present one provides somewhat more information due to the established connection with 
combinatorics and geometry. Also, I was informed about a recent work by Pule and Zagrebnov |9j in 
which inequality (i) is derived and used. 
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